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Abstract 

We propose an accurate method for counting states of close to zero and positive 
cosmological constant in the Bousso-Polchinski Landscape. This method is based on 
simple geometrical considerations on the high-dimensional lattice of quantized fluxes 
and on a probabilistic model (the "random hyperplane" model) which provides a dis- 
tribution of the values of the cosmological constant. Justification of the assumptions 
made in this model are offered by means of numerical experiments. 

1 Introduction 

One of the recent proposals to solve the cosmological constant problem in cosmology is pro- 
vided by string theory. By dimensional reduction from M-theory to 3+1 dimensions, vacua 
of the effective theory are classified by means of a big number of quantized fluxes leading to 
an enormous amount of metastable vacua, the Bousso-Polchinski (BP) Landscape [lj. The 
cosmological constant problem, namely the smallness of the observed vacuum energy den- 
sity in the universe [2j [3] , can be solved by the presence in this model of a huge number of 
states of very small, positive cosmological constant, together with a dynamical mechanism 
given by eternal inflation [I] which allows to visit all the vacua. An anthropic selection is 
then advocated to explain the smallness of the observed cosmological constant (5j E] . 

In order to quantify this selection a counting of accesible states in the Landscape is 
needed. Two ways of counting have been introduced so far: 

• The simplest Bousso-Polchinski count, which computes the volume of a spherical 
shell of small thickness in flux space and divides it by the volume of a cell. 
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• The entropy count of Bousso-Yang, which computes the entropy of the occupation 
number of each flux assuming that they are independent. 

In the following subsection we will briefly review these two counting methods, and drawing 
on them we will propose an alternative one. 



1.1 The Bousso-Polchinski count 



The first estimate of this number is given by the BP count. We will now review this 
argument (see [T]). A vacuum of the BP Landscape is a node in a J-dimensional lattice 
C generated by J charges q%, ■ ■ ■ ,qj determined by the sizes of the three-cycles in the 
compactification manifold. The lattice C is 



£ = {(niQi, ■ ■ ■ ,njqj) G 



Til, 



rij e z} . 



(1) 



The j-th coordinate of a point in the lattice is an integer multiple of the charge qj, and 
therefore a vacuum is characterized by the integer J-tuple n = (m, • • • , nj). 

A fundamental cell (also called Voronoi ce/Q) Q n around a node n in a lattice C is 
the subset of 1R J which contains the points which are closer to n than to any other node 
of C. Thanks to the discrete translational symmetry of our lattice ([I]), all fundamental 
cells in £ are translates of the fundamental cell around the origin Qo = Q, which we can 
parametrize in Cartesian coordinates as a product of symmetric intervals 



x eQ 



x 



X J, 



with Xj G 



2 ' 2 



i.e. 



(2) 



(3) 



The cosmological constant of vacuum n in the BP model i^] 

1 J 

K{n)=K + -Y,n)q 2 r 



(4) 



In Q, A is an a priori cosmological constant or order —1. Each value of A > A defines 
a spherical ball on the J-dimensional flux space of radius Ra = a/2(A — Ao). We call this 
ball £> J (A). We take small values of the charges qj (natural values expected by BP are of 
order |) in such a way that the ball can contain a huge number of fundamental cells. 

The number of states in the Weinberg Window, that is the range of values of the 
cosmological constant allowing the formation of structures (like galaxies) needed for the 
formation of life as we know it [5], is obtained by computing the volume of a thin spherical 



1 Also called Wigner-Seitz cell in solid state physics, the Voronoi cell of a point P in a discrete set S of 
a metric space M is the set of points of M which are closer to P than to any other point of S. 
2 We use reduced Planck units in which 8irG = h = c= 1 . 
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shell in flux space (the realization in the BP Landscape of the Weinberg Window) divided 
by the volume of a cell in the lattice: 



WW 
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(5) 
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where R n = R 



A=0 



2 Aq , and the volume of the J — 1 dimensional sphere is 



vol S 



3-1 



2ni 



(6) 



This method can be naively expected to yield a good estimate when the linear dimensions 
of the cell are small when compared to the thickness of the shell; but this condition is not 
satisfied in the BP Landscape. Nevertheless, the result of this counting formula is very 
good when compared with actual numerical experiments. We will derive the true condition 
of validity of the BP count in our own framework in section [2] below. 



1.2 The Bousso-Yang count 

The second estimate is given by the Monte Carlo numerical simulation by Bousso-Yang 
(BY) in ref. [7j. They compute the number of states by means of the Shannon entropy 
of the occupation number distribution of each flux in a sample of states obtained in two 
ways: 

• Drawing the states from a canonical ensemble distribution with spherical symmetry. 

• Drawing the states as the output of a decay chain using a dynamical relaxation 
mechanism. 

This method has an advantage: it incorporates the dynamical relaxation mechanism of 
Brown- Teitelboim (BT) [HI E], thus quantifying the dynamical selection effect, but two 
drawbacks should be mentioned: 

• It makes it necessary to choose particular values for the charges, thereby providing 
no explicit dependence of the computed number of states with the charges or the 
dimension, and 

• They assume that the states of low positive cosmological constant (referred to as 
penultimate states) have fluxes which are considered independent random variables. 
A hypothesis which they accept as not true but which they suggest how to correct. 

We will briefly review their use of the Shannon entropy for counting. 
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For any subset E having A/"e nodes in the Landscape (E can be the set of penultimate 
states, for example), we can define the uniform probability over E as follows: 



/const if n e E, 
P E (n = n = <^ 7 
10 if n f E, 

where n is a random variable which can take values over the whole Landscape C viewed as 
a subset of Z J , that is, n can be any integer J-tuple n = (n 1; ■ ■ ■ , nj) e£cZ J with equal 
probability, namely The Shannon entropy S E of the uniform distribution P E satisfies 

A/s = e 5E = exp ^ P s (n = n) log P s (n = n) , (8) 



nez 



as can be seen by subtituting (17) into (8) taking const = jj^- If Ps were not constant over 
E, then it must have support larger than E in order to satisfy eq. ([8]). 
If the fluxes were independent, this joint probability would split: 



j 

Pdn = n) = l[P j (n j = n j ), (9) 

3=1 

and therefore the correspondent Shannon entropy would be simplified to 

4 ndep = - EE p ffi = n ) log p ^ = n ) • ^ 10 ) 

j=l neZ 

Unlike P^, the distributions Pj{nj = n) are much simpler to estimate by sampling a small 
portion of E. But, as long as they are not constant, its support covers a much larger 
region than E, and furthermore this region has the symmetry of the cell Q rather than 
the symmetry of E, so we can expect the simplified entropy S^ dep of eq. (10) to be much 
larger than the true entropy S^. 

A numerical experiment can be illustrative. Needless to say, the only way to compute 
the correct uniform probability Ps is to exhaustively compile all the elements of the set E, 
and this is not possible over a realistic Landscape. So we have taken a very simple model 
of J = 3 fluxes with A = -1 and charges q 1 = 0.02988, q 2 = 0.04988, g 3 = 0.06988, and as 
the set E we have considered the secant states (see next subsection). A brute-force count 
of this set is easily carried out by enumerating all states in it, resulting in S% = 9.79222 or 

A/e = 17,894 states. (11) 

The same enumeration of all states allows us to compute the three probability distributions 
Pj{nj = nj), and from these we obtain the additive entropy S^ dep = 12.27932, and a state 
count of 
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which, indeed, means an enormous difference. 

It must be stressed that the goal of the BY work is to quantify the dynamical selection 
effect of the BT decay chain, and they do this by computing the quotient of the number of 
states with and without dynamics. While both numbers are affected by the independency 
error, their quotient may be free of errors. We are interested in studying this point in our 
future research. 

1.3 Our counting method 

Our proposal is based on the following kinds of states one may encounter near the null 
cosmological constant surface in flux space: 

• Boundary (or penultimate after BY) are those states in which a BT decay chain 
can end before jumping into the negative cosmological constant sea. So we define a 
boundary state as one having 

(1) positive cosmological constant, and 

(2) at least one neighbor of negative cosmological constant. 

• Secant states have the property that their Voronoi cell in flux space has non-empty 
intersection with the null cosmological constant surface in flux space. Note that a 
secant state may have negative cosmological constant. 

These two categories are not equivalent; a boundary state may not be secant if it is far 
enough from the null cosmological constant surface, and a secant state may not be boundary 
if it has negative cosmological constant. So we are interested mainly in the states which 
are both secant and boundary, because all the states in the Weinberg Window are in this 
category. Figure [T] illustrates the differences between these two kinds of cell. 

Our strategy would be as follows. We will count the states in the Weinberg Window 
using the following elementary formula: 

A/ww = \Ar s P(A G WW) , (13) 

where 

• A/ww is the number of states in the Weinberg Window, 

• As is the total number of secant states, 

• the | factor is the (first-order) approximate fraction of positive cosmological constant 
secant states, 

• P(A G WW) is the probability that a random secant state has a positive cosmological 
constant in the Weinberg Window. If we call a number of the size of 10~ 120 Aww, 
then 

P(A G WW) = P(0 < A < A WW ) , (14) 
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Figure 1: States are shown on the A = surface in a J = 2 BP Landscape. We see secant 
states which are not boundary and viceversa (a), regions in which S:B are in a 1:1 ratio 
(b) and regions in which S:B are in a 2:1 ratio (c). 



where the probability must be computed using the distribution of the cosmological 
constant as a random variable over all the secant states. 

In section [2] we first compute the number As and check it with simple models with two or 
three fluxes, where brute-force counting is feasible. In the following section we propose a 
probabilistic model (and check it against numerical data) which leads to the distribution of 
the values of the cosmological constant restricted to the secant states (section |3]). We use 



this distribution to compute the probability (14) in sectionHl Surprisingly, our results show 



that this more precise counting method yields the same result as the BP count. Finally, 
we summarize the conclusions in section HI 



2 Counting secant states 

We start with the observation that the number of secant states which we can find by looking 
from the origin in a given direction is not a constant, as can be seen in the histogram shown 
in figure [2j In this figure, we can also see the theoretical density of states which we now 
derive for the J = 2 case. 

Let N(6) be the number of secant states on the first quadrant of a BP Landscape 
which are between the 1-axis and a straight line drawn at an angle 9, as in figure [3| The 
number of states along the arc between the 1-axis and the #-angle line can be accurately 
approximated by the length of the segments which the states along the A = circle are 
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Directional density of secant states (J = 2) 




Figure 2: Directional density of states in the BP Landscape shown in figure [I] computed 
numerically compared to the continuous version v{6) of equation (16). All plots in this 
paper regarding statistical analysis were done using R [TO] . 




7 



covering. These lengths can be straightened as shown in figure |3j to yield (here we call 
R = Rq = v/2|Ao| for convenience) 

R^0 + R(l-cose) 
12 qi 

Its derivative is the directional density of states: 

diV /COS6 1 sin^x R , „ , Rq ■ v . . 



In equation (16), volQ = ?i?2, ^ = (cos^, sin^) and q = (gi, q^)- Only the first quadrant is 
considered here and the absolute value on the components of v must be taken if we want 
to extend equation (16) to all quadrants. The formula thus obtained is plotted in figure 
[2] over the numerical data, giving an accurately approximated density. Of course we can 
estimate the total number of secant states by 

A/l hco = 4 P i/(0) d6 = AR{qi + g2) . (17) 
5 Jo Q1Q2 

Using the J = 2 values q x = 0.01494, q 2 = 0.02994 and A = -1, we have N^ heo = 
4 x 141.8945 = 567.5779, being the correct number, obtained by an exhaustive search, of 
143 states, 2 of them on the axes, which gives a total of 

A/£ um = 4 x 141 + 4 = 568 « 567.5779 = A^ heo , (18) 

The two results clearly agree. 

For a model with J = 3, we can see that the Voronoi cells of the secant states on the 
first octant project their faces over the first quadrants of planes 1-2, 1-3 and 2-3. The area 
of these quadrants is |i? 2 and the number of states needed to cover this total area is 

Arthpn irR 2 ttR 2 ttR 2 ttR 2 , N 



Let us now take the density (16) and extend it to J = 3. The correct scaling will be 

(20) 

where v is a norm-one vector on the S 2 sphere. Let us integrate this function on the first 
octant: 

1 f R C~2 

- / v{v) (\.Vt 2 {v) = / sin6>d# / d0(gi sin 6* cos 0+g 2 sin 9 sin 4>+q^ cos 9) . (21) 

8 J& qiqzq?, Jo Jo 

The ^-integrals have the values 1, 1, |, and the ^-integrals have the values |, |, |, so that 
all integrals coincide and the result is (19). 
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Using the values q l = 0.01494, q 2 = 0.02494, q 3 = 0.03494, A = -1 for the J 
model, we have an 8-octant total number of secant states of 



JQ- = ^- (qi+m + q3 ) = 4 * x °f 482 5 = 72220.02 . (22) 
5 4volQ vyi y H3J 1.301877- 10- 5 v ; 

A brute- force search of the secant states gives 9205 states on the first octant. But 3 are on 
the axes, 96 + 135 + 152 = 383 on the coordinate planes and the rest (9205 — 383 — 6 = 8819) 
outside, so the total number of secant states found in this model is 



7\/^ um = 2x3 + 4x 383 + 8x 8819 = 72090 . (23) 

The agreement in this case is not complete but it is satisfactory. 

The discrepancy may be caused for the sampling method of the secant states in the 
brute-force search. It consists on uniformly sampling the Sq sphere of radius a/2|Ao| 
correspondent to A = 0; each point thus sampled hits on the intersection between the 
sphere Sq and the Voronoi cell Q n of certain secant state n. Therefore, the probability for 
this secant state n to be selected by this method is proportional to the area (( J— 1)- volume) 
of the intersection: 

p ( „H«> (24) 

If we call the set of secant states S, the complete and disjoint partition (tessellation) which 
the Voronoi cells induce on the sphere guarantees that this probability is normalized: 

5>(n) = l. (25) 

Thus, only the states with bigger intersection area will be selected in a sample taken by 
this method. In addition, the intersection volume is positive but it does not have a positive 
lower bound, so in principle there can be arbitrarily small intersections which will not be 
detected by this method. In the J = 2 case, 10 4 points were needed to find the existing 143 
secant states (see figure JTJ, but using 10 8 points on the J = 3 case yielded 9134 states of 
the 9205 revealed in the 10 9 sample. If the set S contains states with a probability below 
10~ 9 , these will not be found using this method. So the agreement may increase by taking 
a sample size greater than 10 9 . 



Let us extend the directional density of states (20) to an arbitrary dimension J and to 
all "quadrants": 

^ ) = voTQ 9>1 ' (26) 

where v = (v\, ■ ■ ■ ,vj) G S 1,7-1 , \v\ = (\vi\, • • • , \vj\) and q = (qi, ■ ■ ■ ,qj). The number of 
secant states is 

N s = f u(v)dn J _ l (v) = I ^—y2q i [ \vi\dnj-iiv). (27) 

JSJ-l VOl Q ^ Jgj-l 
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The integration measure on the sphere dQj-i(v) is invariant under SO (J) rotations, and 
there is a rotation which can transform a given coordinate V{ in the simplest of them, vj, 
in the following choice of coordinates for the sphere: 



v = w sin 9 + e j cos 9 , w G S J 2 , 9 G [0, ir) , 
where ej is the last vector in a R J basis. The integration measure reduces to 



dQ 



sin J " 2 9 d9 dtt 



J-2{W) , 



so that 



\vi\ dTLj-Av) 



i s j-i 
We have then 



cos9sm J ~ 2 9d9 



dQj- 2 (w 



2 vol S 



J-2 



S J-2 



J -I 



Ms 



2R J - 1 vo\S J ~ 2 



j 



2 J R J " 1 gvolS 



J-2 



(J-l)volQ j^ 1 ' 1 voig 



(2J 



(29) 



(30) 



(31) 



In a somewhat nonstandard way, we have defined the quantity 

j 

i i 



1 J 

i=i 



ft = ?■ 



(32) 



J-2 



Equations (|22|) and (|17|) are special cases of (|31|). A more explicit form of (|31|), using 

2 



i? = a/2|A | and vol S 



27r^/r(^) is 



nil ft 



(33) 



Let us now derive the condition of validity of formulae ( |31| |33J) . Note that we are summing 
the number of times that a fundamental cell fits in a ( J — l)-quadrant; in order for this 
number to represent the actual number of secant states, the excess ( J — 1)- volume of the 
cells trying to fit the boundary must be very small when compared with the (J — 1)- 
volume of the ( J — l)-quadrant. The ( J — l)-volume of the cell projection is different in 
each direction, so we can compute a kind of mean value of the projection volume as 



j-i 



j j 



J 



j = l 8=1 

#3 



q J 



(34) 



where q and qn are respectively the geometric and harmonic means of the charges. This 
mean projection volume must be much less than the (J — l)-volume of a single (J — 1)- 
quadrant, that is, 

volS J - 2 . (35) 



« 



2- 7 - 1 J - 1 
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After substituting R = y/2\A Q \ and vo\S J - 2 = 2vrV/r(^), we fincj^ 

u 2 

Jy— < vre « 8.539734 . (36) 

3 The random hyperplane model 

We assume two basic features of the secant states in a BP Landscape: 

• The Voronoi cells of the secant states are small enough to replace the A = sphere 
which intersects the cell by its tangent hyperplane. 

• The orientations of the hyperplane intersecting the cell are random if one picks a 
secant state at random. 

In this way, we propose to study the set of secant states by choosing a probability measure 
on the secant hyperplane set. In this section we first parametrize the hyperplane space, 
and then we choose a probability measure on it. 

3.1 The hyperplane space 

We define H as the set of all hyperplanes in IR J and Hq as the set of all hyperplanes with 
non-empty intersection with Q: 

H Q = {h G H: Qf]h^ 0} . (37) 

In order to specify a hyperplane h G Hq we must choose, first of all, a normal vector of 
unit norm, v G S 7 ^ 1 . A point in h must also be given; but there is an infinity of possible 
choices here, and we must provide a unique prescription, for example the point z G h 
closest to the origin O (chosen here as the center of Q). Note that this point z may lie 
inside or outside Q. 

The pair (z, v) has a lot of redundant information, because the directions of the vector 
Oz = (zi,-- - ,zj) and v = (vi, ■ ■ ■ ,vj) coincide, that is, Oz = ±\Oz\v. Therefore, we 
must retain only the norm of Oz and the whole v; letting p = \Oz\, we can identify h with 
the pair (p,v). Also, we must take into account the fact that v and — v represent different 
hyperplanes if p ^ 0, but the same hyperplane if p = 0, so the vector v is defined up to a 
sign in the case p = 0. 

An alternative equivalent prescription is to consider directions up to a sign for all p 
(the point v is on S' 7 ' 1 with antipodal points identified, i.e. a point on the projective space 
P J_1 ). In this case, p represents not only the minimum distance from the hyperplane to 

incidentally, this relation resembles the so-called t'Hooft limit in field theory, in which the number N 
characterizing the gauge group tends to infinity and the Yang-Mills coupling constant <7ym vanishes with 
the product vV<?ym (^ ne t'Hooft coupling) held fixed. 
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the origin O, but also the position of the hyperplane "above" or "below" the origin. We 
find the former point of view more adequate. 

The domain of definition of p is a direction-dependent positive interval I v = [0,a(v)], 
and is determined by the restriction that the hyperplane it represents has non-empty 
intersection with Q. In particular, the point z of minimum distance can lie outside Q; 
thus, the "hyperplane space", when compared with "physical" space, comprises a larger 
region than Q. 

We now compute cr(v), defined as the maximum distance to the origin of the closest 
point to the origin of a hyperplane with non-empty intersection with the cell Q. Clearly, 
the corners of the cell are the most far away points in Q, so for each direction v, the last 
(most distant) hyperplane orthogonal to v and with non-empty intersection with Q must 
contain one of the 2 J corners of the cell. The equation of such a hyperplane is 

v{x-c v ) = 0, (38) 

where c v = |(si<7i, ■ ■ ■ , sjqj) is the unique corner out of 2 J which belongs to the same J- 
quadrant as v, and the Sj are signs ±1, indeed the same signs of the components of v, that 
is, v = (si\vi\,--- ,sj\vj\). This hyperplane is already in normal form, so its minimum 
distance to the origin is 



cri v) = v ■ a 



2 S ?* s i IV* I = 2 S q i\ v i \ = \l-\ v \- ( 39 ) 

3=1 3=1 

Given a direction v, the closest point can be found as the intersection of the line x = vt 



with the hyperplane (38). This method is used to construct the hyperplane space for the 
J = 2 case in figure |4j where it is compared with formula (39). The function cr(v) can 



be plotted for J = 2 and q\ = 0.01494, q 2 = 0.02994, versus a polar angle 9 defined as 
v = (cos 9, sin 9), see figure |5j 



3.2 A probability measure in the hyperplane space 

Now we need to define a probability measure on the hyperplane space just parametrized. 
The simplest choice is the uniform probability on Hq\ 

dP(h)= dP(p,v) = K- 1 x HQ (h)dpdQj^(v), (40) 

where XH Q {h) is the characteristic function on Hq, whose effect is simply to restrict the 
integration domain to Hq, and dQj^i(v) is the volume measure in S J ~ l . K is a normal- 
ization constant: 

K= I dpdQj^(v). (41) 



We can compute K by first integrating out the p variable in (41): 

f(T(l>) 



K= I [ dpdn J „ 1 (v)= [ <7(u)dfij_i(u). (42) 
Js 1 - 1 Jo Js'- 1 
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Figure 4: Construction of the hyperplane space in the J = 2 case. The blue contour 
comprises all points which represent a different hyperplane intersecting the cell. Here, the 
corner A is used to construct the first quadrant part of the contour, and a = \OB\. The 
green points have been generated using formula (39) on the third quadrant. 




Figure 5: The function a in the J = 2 case, with q\ = 0.01494, q2 - 
(cos 9, sin 9). The cell boundary is also displayed. The corner is at 9 Q 
the maximum of a. 



0.02994 and v = 
- arctan — , and is 

91 ' 
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By comparing equations (26) and (39) we observe the following relation between the max- 
imum distance of a hyperplane in a given direction o~(v) and the directional density of 
states v{v): 

. . 1 vol Q . . , . 

<v) = --^u(v), (43) 

so that using the definition of the density of states 

M s = [ v(v)ttlj- x {v) (44) 
Js 1 - 1 

we have 

K= f a(v)dQ J ^(v) = A ^^ = qvo\S J - 2 . (45) 



In the last formula we have substituted the expression for Ms found in equation (31). 

The main reason for choosing the uniform probability measure is its simplicity. The 
physical reason is that among all measures in a given compact space, the uniform one 
has maximum Shannon entropy. In the following, we will justify our choice by more 
quantitative means. 

First of all, the marginal probability in the v variable associated to the uniform measure 
is proportional to the directional density of states v(v), with a normalization-to-one factor: 

f f 7 ^ 2E J ~ l i/(v) 

/ dP(p, v) = K- 1 dQj^(v) / dp = -= — a(v) dQj^(v) = dSlj^v) . 

Jpe[o,(T(v)] Jo JMsVOiLJ J\ s 

(46) 

Therefore, the uniform probability reproduces the correct directional density of states as 
observed in the numerical experiments. 

On the other hand, we can plot instances of the set of secant hyperplanes and compare 
them with simulated uniform points in hyperplane space. The simplest of such plots is 
for J = 2; with A = -1 and charges q x = 0.001494, q 2 = 0.002994, we show the (9, p) 
points from the actual secant hyperplanes versus a simulated set of uniform points of equal 
sample size in figure [6} 

The geometrical nature of the secant hyperplane sample is revealed in the structures 
shown in fig. [6] (left), which introduces correlations in the spatial sequence. In contrast, 
a uniform sample in hyperplane space (fig. [6j right) shows no correlations, and the only 
structures which we can see are voids and clusters bigger than the ones present in the 
former case. 

Both samples thus obtained (the secant hyperplane set, or geometric sample, and the 
uniformly distributed points, or simulated sample) cover the hyperplane space, so we can 
approximate the distribution of the former set of points by a uniform probability, thus 
neglecting the spatial correlations. 

But the efficiency of the covering is not equal in both cases: the geometrical sample 
shows smaller voids and clusters than the simulated one. This can be seen by counting the 
number of points inside a circle of random position and given radius inside the hyperplane 
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Secant A > hyperplanes in a J = 2 Landscape 



Uniformly distributed points in J = 2 hyperplane space 




Figure 6: Actual secant hyperplane plot of a J = 2 Landscape in hyperplane space (ge- 
ometric sample, left) and a uniformly distributed sample in the same domain (simulated 
sample, right). Note the structures in the former, and the bigger voids and clusters in the 
latter. 



Secant A > hyperplanes in a J = 2 Landscape 




8/1 

700 points, q, = 0.001 494, q 2 = 0.002994 



Figure 7: Illustration of the circle method to measure the voids and clusters distribution. 
Small circles are thrown randomly inside hyperplane space. Green points are the sample 
points inside the circles; red points are the random centers of the circles. 
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Histogram of number of points (secant sample) 



Histogram of number of points (uniform sample) 




5 10 15 5 10 15 

Number of points inside a circle of radius 0.05 Number of points inside a circle of radius 0.05 



Figure 8: Histograms of the number of points inside a random circle in the geometric (left) 
and simulated (right) cases. Both are well described by a Gaussian; for a comparison, see 
text. 



space in both cases (fig. [7] illustrates this "circle method"). When this is done for a 
big number of such circles, the voids and clusters induce a fluctuation in the number of 
"inner" points, which will be greater if the inhomogeneities are bigger. Fig. [8] shows the 
result presented as two histograms of the number of inner points for 10 3 circles. These 
can be well approximated by Gaussians, and we can see that the variance of the simulated 
sample is bigger than the variance of the geometric one: that is, the geometric points 
cover hyperplane space with a more regular (less random) pattern, in such a way that each 
point tries to avoid the clusters and tends to fill the gaps, a characteristic behavior of the 
so-called quasi Monte Carlo (qMC) sequences [TT] . 

Putting the comparison in more a quantitative fashion, the mean and standard devia- 
tion of the geometric and simulated samples are 

/i Geo = 6.69697 , /i Sim = 5.787115 , 

a G co = 2.185198 , a Sim = 2.474453 . ^ ' 

The different mean values with the same number of points indicate that the density of 
points is not the same in both samples. This happens because the spatial correlations 
in the geometric case make the left and right sides underpopulated (see fig. [6]), so that 
the effective volume covered in this case is smaller. We can take this into account by 
considering the quotient between the standard deviation and the mean as a measure of the 
clustering of the data, which is an effective scaling of the data to have mean value 1. This 
leads to 

a ^ = 0.3262966 , — = 0.4275798 . (48) 

/^Geo /^Sim 

This means that the fluctuation in the point density is greater for the simulated sample. 
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The relative difference is considerable ( a4275 7 ^~° 9 g 62966 = 0.3 104 23, that is, 31%). This 
indicates that the covering of the hyperplane space by the geometric sample is much more 
efficient than the covering made by the simulated one, and thus its randomness is much 
smaller. 

Summarizing, the approximation made by assuming uniformity in the distribution of 
secant hyperplanes is equivalent to neglecting the spatial correlations present in the ge- 
ometric case, and we can expect this approximation to work well because of the more 
efficient covering of hyperplane space in this case. 



4 Number of states in the Weinberg Window 

The next step in our calculation is to compute the distribution of values of the cosmological 
constant. The random hyperplane model allows us to do this by means of the marginal 
probability distribution in the p variable, whose density will be called ui(p): 



f 



ves- 1 - 1 



dP(p,v) = u(p) dp. 



(49) 



Using the following relation between p (the minimum distance from a tangent hyperplane 
to their secant state), R = a/2|A | (the radius of the A = sphere) and the Euclidean norm 
of the secant state in flux space (which in turn is related to the cosmological constant) 



R + P 



A 



1 J 



2 2 



t=l 



p(A) = v/2(A-Ao)-v/2|A c 



(50) 



the A distribution can be found once the function cu(p) is computed: 

/(A) dA = ou[p(A)]^- dA = u[p(A)] 

JK ' [PV ;j dA lPV n ^2(A - A ) 



(51) 



The u> function is easily computed once the integration of the characteristic function in 



the (p, v) variables is reversed. The p variable takes a maximum value <r max = |y Yli=i ih 
and the v variable can range across 5 J_1 . So we define the integral of an arbitrary test 
function <p(h) on the hyperplane space Hq against the probability measure dP(h) as 

/ XH Q (p,v)<f>(p,v)dpdQj-i(v) . 

'(p.^GlO^maxlxS- 7 - 1 



(f)(h) dP(h) 



K 



(52) 



For the iterated integrals, we have the following identitjj^] (Fubini's theorem): 



dp / dVtj-i{v) (p(p,v) 



s J - 



dfij_i(t;) / dp <f>(p,v) 



(53) 



In equation ( 53 1 , we drop the constant K computed in ( 45 1 . 
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where the sets I v C [0, cr max ] and J p C S J 1 appear owing to Xh q , and are defined as 

I v = { P E [0, ow] : (p, v) E Hq] = {p E [0, a max ] : < p < a{v)} = [0, a{v)\ , 
Jp = [v E S J - 1 : (p, v) E H Q ) = {v E S J ~ l : a{v) > p) . 

Note that the restrictions in both sets are the same, p < cr(v), but in X v the direction v 
remains fixed, and in J p the distance p remains fixed. So, X v is simply the interval [0, o[v)\ 



used in the integral (46), and J p is a set which is the whole sphere when p < | min{gj} 
and a 2 J -connected-component set when p > \ maxjg,'}. If we carry out the v integration 
of dP(p, v) but we stop before the p integration, we have 



■' Jp 



To compute the integral in (55) in closed form is not an easy task, and even if it were, 
it would be useless, because the function u is built by gluing polynomial functions in a 
J-dependent number of intervals. 

Let us compute uo for the simplest J = 2 case. The set J p is here an interval, which 
can be computed from the a function which represents the upper boundary of hyperplane 
space, see fig. [5] A horizontal line at height p has two, one or no intersection with o{v) 
in the first quadrant depending on p. Taking qi < q 2 , there are no intersection points if 
p < y, there is a single intersection point if y < p < y and there are two intersection 
points if y < p < cr max . The equation to be solved for the intersection points in the first 
quadrant is (remember that in J = 2 we have v = (cos 8, sin 8)) 

p = ^ [qi cos 8 + q 2 sin#] . (56) 



Setting x = cos8, we have a quadratic equation whose solution has two branches: 

x± 



a 2 

max 



2 P± ■yV / ° r max-P 2 J • ( 57 ) 



These branches only make sense for p > y (for x + ) and for p > y (for x_). We define 
$-1- = cos -1 x±, so we have, in the first quadrant only, 




The uo function is, then, 




ifp<f, 

if f < p < f , {51 

if f < P < | + ?2 = °"max- 



ifp<f, 

if f < p < f , (59) 

if f < p < lVtf + ti = 



a 



max- 
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Histogram of the p parameter in a J = 2 Landscape Exact and numerical evaluations of the co(p) function (J = 2) 




2 Landscape. The solid line is the uj 
the points represent the Monte Carlo computation of 



Figure 9: Histogram of the p parameter in a J 
function computed using eq. (59), 
the uj integral (left). Relative error between the numerical and exact evaluations of the uj 
function in the J = 2 case. The MC integration used 10 7 points in the circle S 1 (right). 



With J = 2,K = 2(q l +q 2 ), and taking A = -1 and charges q x = 0.001494, q 2 = 0.002994, 
we can plot this function and compare it with the true p histogram of the secant states; this 
is done in figure [9] (left), where we see big fluctuations of the number of states versus the 
theoretical uj function. This strong fluctuation is a consequence of the relative sparsity of 
this Landscape. The small sample size makes the histogram oscillate around the theoretical 
curve; we believe that this effect is only present in J = 2, as can be seen in the subsequent 
figures. 

In figure [9] (right) we can see the relative error between the exact formula (59) and the 



Monte Carlo evaluation explained below (eq. (61 )). They agree in the constant region, and 
the error remains smaller than 10~ 3 for all values. 

We can compute the values of the uj function numerically by rewriting it in the following 
way, using the unit step function 9(x) for restricting the integrand to J p and the formula 



(45) for the normalization constant K: 



vol^- 1 f 1 diWfo) 

= q^oiS^ J SJ -, 6[(J{V) ~ p] ^oLS^ • ( 60 ) 

Thus, the formula splits in a constant times the mean value of the step function with respect 
to the uniform probability measure on the sphere. This mean value is easily computed using 
a simple Monte Carlo technique, that is, sampling the unit sphere with a large number N 
of points t>W, we have 

vol S^ 1 1 - 
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The Monte Carlo evaluation is well suited for this task, because the integrand is bounded 
and the integration domain is compact. 



In fig. 10 we can see the MC-computed uo distribution for Landscapes with J = 3 and 
J = 4, where brute-force data are available. The J = 3 oscillates much more than the 
J = 4 one, due in part to the bin width, which contributes to make the histogram smoother 
with a bigger sample size. But the mean value in the constant region (as well as the global 
fit to the entire theoretical curve) seems to be better adjusted in the J = 3 case. The reason 
for this apparent disagreement between the brute-force data in J = 4 and the theoretical 
MC-computed to curve is the sampling method used to find the secant states in this case. 
This sampling method, described above (see section [2]), has the disadvantage of missing 
the states with very small area of the intersection polytope between the A = sphere and 



the Voronoi cell of the secant state. Using eq. (31 ), we can estimate the number of missing 



states; with A = -1, q x = 0.01494, q 2 = 0.02244, q 3 = 0.02994, g 4 = 0.03744, we have 

A/"j h j° 4 = 6, 605, 383 secant states. (62) 

The brute force calculation yielded 406,715 states in the first 4-quadrant; taking into 
account degeneracies (the degeneracy of each state is g = 2 4 ~ z , being z G {0, 1,2,3} the 
number of zero-flux components), we obtain 

^brut^force = 6; 2 45, 948 secant states. (63) 
Their difference relative to the mean degeneracy is 

^/y"thco ^y*brute-force 

— 5 ' J ~ 4 — l5 ' J ~ 4 = 23,405.19 first 4-quadrant secant states. (64) 

(9) 

These states amount to 5.4% of the total, and they lie completely in the high-p region. 
This is the reason why the theoretical curve disagrees with the histogram in the J = 4 case: 
the histogram is normalized to have total area one and it misses 5% of states of higher p, 
so we must conclude that the data range of the histogram must be shorter and its constant 



region must be higher, as shown in fig. 10 



But we don't need the numerical evaluation to analyze the to function. From eq. (59) 
we can see that when p < | min{gj}, the value of the step function is always 1, so we have 
the exact result 



volS 7 - 1 v^Fr(^) r 1 r 1 

U{P) = q^oiS^ = qT(i) fOT P<2 mm ^}' ^ 

that is, the distribution of p values is exactly constant at small p, p < | min{gj}. Beyond 
this value, the step function finds regions where &(v) < p, and the mean value begins to 
decrease monotonically. When p > cr max , no point in the sphere is captured by the step 
function, and u vanishes. We have checked this general behavior on high-dimensional BP 
models before and after performing the BT decay chain and this qualitative profile is quite 
robust EL2]. 
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Histogram of the p parameter in a J = 3 Landscape 



Histogram of the p parameter in a J = 4 Landscape 




Figure 10: Histograms of brute-force search for secant states compared with the MC- 
computed uj curve for Landscapes with J = 3 (left) and J = 4 (right). See text for an 
explanation of the differences. 



As the authors of [32] point out, actual histograms of A on particular BP Landscapes 
show a "staggered" behavior which gets smoothed when the bin width increases, but the 
mean value of these oscillations is given by uj. Note that A is really a discrete variable, and 
uj is the density of a continuous one. Furthermore, note that in [13J the authors compile 
statistics on entire Landscape instances, while our statistics refer to the positive A secant 
state sector only. 

So we find that uj behaves as a Fermi-like distribution on the p values, with a medium 
filling level and a decreasing interval which depend on the values of the charges. Further 
quantitative analysis must be done trough the numerical estimation procedure described 
above. 

Now, we can compute the probability for the cosmological constant to lie in the Wein- 



berg Window, using the properties of uj and formula (51): 



P(0 < A < A 



WW 



Aww 



/(A)dA 



volS 



J -I 



Aww 



dA 



qvo\S J - 2 J v^A-Aq) 



vol S J - X A 



WW 



R q vol S 3 2 



(66) 



After inserting this in formula (13), together with eq. (31), we have 
1 2R J - 1 qvo\S 



:J-2 



WW 



X 



volQ 



vol S ,J-1 A W w 
Rq vol S J ~ 2 



R J ~ 2 vo\ S^A 



WW 



volQ 



(67) 
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Number of states in a shell vs. BP count 




^ 1 1 1 1 1 1 — 

0.00 0.01 0.02 0.03 0.04 0.05 0.06 

Shell width (100 values in range) 
q, = 0.02988, q 2 = 0.04988, q 3 = 0.06988 



Figure 11: Comparison between an actual count of the number of states in a shell of 
variable width (stair-like line) and the number predicted by the BP count (plain line). 



This result exactly coincides with the BP count (eq. ([5]) 

One may wonder if the approximations made by the random hyperplane model are 
too crude to achieve a better result than the BP count, or if the BP count is a better 
result than expected. In order to clarify this point, we have again performed a numerical 
experiment to compare the number of states predicted by the BP (and our) count with a 
brute-force search of states in a shell of variable width (that is, varying A-yvw) in a J = 3 
Landscape with 2,333 secant states in the first octant. In our opinion, the nature of the 
example chosen would render the approximations made too crude, but we find a remarkable 
agreement, as shown in fig. 11, and in a range of shell widths much wider than expected. 

We are forced to conclude that the BP count has proved to be a succesful way of counting 
states in the BP Landscape, in fact better than expected. The geometric interpretation 
is curious; if we sum up the volume of the Voronoi cells of the secant states whose center 
are inside a shell of prescribed width, we obtain the volume of the shell with very good 
approximation. In figure 12 we illustrate the spatial location of the states inside the shell 
for chosen shell widths in the J = 3 Landscape of the previous example. These locations 
are uncorrelated for low values of the shell width. 
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Figure 12: Blue lines represent the tessellation induced by the Voronoi cells of the secant 
states in the first octant (whose boundary is depicted by a black line) of the A = sphere. 
Each tessella corresponds to a secant state. The red tessellas have centers inside a shell 
of width 0.001 (top left), 0.003 (top right), 0.01 (bottom left) and 0.07 —all secant states 
with A > — (bottom right). 
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5 Conclusions and future directions 



We have modeled the BP Landscape by means of the random hyperplane model (RHM), 
validated by several numerical experiments like brute-force search and precise state count- 
ing. The RHM provides us with a distribution of values of the cosmological constant and 
reproduces the BP count in a very different way. The numerical experiments performed 
suggest that the BP count is much more precise than expected, far better than the BY 
count, for example, and applicable to a wide range of parameters of the Landscape. 

Further corrections to this formula can take into account the asymmetry of the A > 
and A < states due to the curvature of the A = sphere (neglected here), which is 
a second order effect which depends on the precise values for the charges. For us, it is 
certainly more difficult to take into account the spatial correlations, which would require 
a completely different approach. 

All these considerations lie in the context of finding an a priori distribution of the 
cosmological constant, without any reference to the dynamics of the BT decay chain. In 
our future research we will try to incorporate the dynamical effects as well, in order to 
count also the states selected by dynamics and to quantify the dynamical selection effect 
within this framework. In this respect, we have checked by extensive numerical simulations 
that the Fermi-like profile of the u ditribution on the secant states persists even on a subset 
selected by the BT decay chain [12]. If this effect reveal a feature of the dynamical selection 
or is caused by the geometric structure of the set of secant states is a point which is worth 
of further investigation. 



A Improving the BP count 



We can slightly improve the result obtained using eq. (13) and the random hyperplane 
model (RHM) described above. The first heuristic approach is to replace the | factor in 



( 13 ) by a measure of the volume of a shell above and below the A = surface in flux space. 



This surface is a sphere of radius R = a/2|A |. We will call the inside region of this sphere 
B J (R) and its volume vol£> J (i?) = V J (R) = ^volS 1,7-1 . We can draw a shell of width e 
above the sphere and another shell of the same width below. Then, the relative volume 
r) of the positive-A region in the shell is given by the quotient of the volume of the upper 
shell and the total volume of the two shells: 

„_ V J (R + e)-V J (R) _ (1 + ^) J -1 (6g) 



v j {r + e) _ v j {r _ e) (1 + ±y _ (1 _ ±y 



In the limit e — > or R — ► oo, rj — > | which is eq. (13), and for the special case e = R 
(when the inner shell completely fills B J (R)) rj = 1 — 2 _J , which is arbitrarily close to 1 
for high enough J. So we are replacing the | by a quantity | < rj < 1, and therefore the 
improvement is small. 

It remains to determine what the relevant width e is. To this end we use the RHM: the 
distribution uj(p) has a Fermi-like profile, and thus we can define a kind of "Fermi level" 
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as the width po of a step function distribution that has the same height ujq = lo(0) as ui 
and is also normalized to one, that is, c^oPo — 1- Thus e = po — ~^ can pl & y the role of the 
effective width of the secant state set. 

Another approach which also takes into account the difference in the volumes of the 



inner and outer shells is to modify the probability measure in hyperplane space (40) and so 
the RHM itself. Note that the p variable is the distance between the tangent hyperplane 
and the secant state, so that R + p is the radial spherical coordinate in flux space of the 
point in the hyperplane closest to the secant state. Therefore the natural modification in 



the probability measure (40) would be 



dP(h) = dP(p,v) = K- 1 XHQ (h)(R + p)- 1 - 1 dpdQj^v) . (69) 

In this alternative viewpoint of the RHM the p variable can be positive or negative, but 
the weights of the two possibilities are different because the marginal distribution in the p 
variable is not symmetric: 



dP(p, v) = K~\R + p) J -^(p) dp , (70) 



where u is the same function as before extended to negative values of p by symmetry 



u)(—p) = oj(p). Now, we change eq. (13) by 



AW = A/" s P(0 < A < Aww) , (71) 

because the new probability distribution carries the difference between the positive and 
negative A states. 

As long as p <C R, this modification will be harmless and the new probability will be 



exactly | the old one, recovering (13). But if p takes on high values, then the asymmetry 



will be enormous: the distribution (70) will develop a sharp peak in the positive p region, 
far away from the Weinberg Window. 

To see how this phenomenon can happen, we can remember that p is bounded by its 
maximum value cr max = ih which behaves as \\J~Jq, being q the square root of the 

second moment of the charges. The extreme case cr max ~ R can be written as 

J^- ~ (2v/2) 2 = 8 . (72) 
l A o| 



But then we are leaving the condition of validity (36), where we have replaced p by q. 

Note that unfamiliar things happen when conditions like (36) are violated. We can 
see this by choosing a fixed small typical charge q <C R and increasing the dimension J. 
The distance from the corner of the cell and its center, \q\/~J ', can thus reach the radius 
R, that is, the Voronoi cell around the origin in flux space will eventually find its corners 
touching the surface of the sphere^] of radius R. At this point all states inside the A = 

5 Thc cell acquires a kind of "dendritic" structure, with its 2J faces lying far away from the sphere 
and its 2 J corners touching it. An enormous fraction of the volume of such a small, J-dimcnsional cell is 
located at its corners. 
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sphere will be secant; of course many states outside the sphere will be secant also, and a 
huge fraction of them will be in the corners of a parallelotope enclosing the sphere. But 
the vast majority of these states will not have neighbors of negative cosmological constant, 
and therefore the typical value of A in these states will be large, hence the peak in the p 
distribution mentioned above. 

We see that once the validity condition is violated, the majority of the secant states 
are no more boundary states, and thus the secant states in this regime have no special 
relation with the Weinberg Window. But the secant state set S is simply a tool to study 
the boundary state set V, which comprises the states selected by dynamics, and we can 
use S instead of V only when they significantly overlap. We believe that this is the case 
when the validity conditions are satisfied. 

Acknowledgments 

We would like to thank Roberto Emparan for useful discussions. We also thank Concha 
Orna for carefully reading this manuscript. This work has been supported by CICYT 
(grant FPA-2006-02315) and DGIID-DGA (grant 2007-E24/2). We thank also the support 
by grant A9335/07 (Fi'sica de alta energia: Particulas, cuerdas y cosmologia). 

References 

[1] R. Bousso and J. Polchinski: Quantization of four-form fluxes and dynamical neutral- 
ization of the cosmological constant. JHEP 06, 006 (2000) [hep-th/0004134] . 

[2] S. Weinberg: The cosmological constant problem. Rev. Mod. Phys. 61, 1-23 (1989). 

[3] R. Bousso: TASI Lectures on the Cosmological Constant. Gen. Rel. Grav. 40 (2008) 
607 [0708.4231 [hep-th]]. 

[4] A. H. Guth: Inflation and Eternal Inflation. Phys. Rept. 333 (2000) 555 and references 
therein [astro-ph/0002156] . 

[5] L. Susskind: The Anthropic Landscape of String Theory, [hep-th/0302219] . 

[6] S. Weinberg: Anthropic Bound on the Cosmological Constant. Phys. Rev. Lett. 59, 
2607-2610 (1987). 

[7] R. Bousso and I. Yang: Landscape predictions from cosmological vacuum selection. 
Phys. Rev. D 75 (2007) 123520 [hep-th/0703206] . 

[8] J. D. Brown and C. Teitelboim: Dynamical neutralization of the cosmological constant. 
Phys. Lett. B195, 177 (1987). 

[9] J. D. Brown and C. Teitelboim: Neutralization of the cosmological constant by mem- 
brane creation. Nucl. Phys. B297, 787-836 (1988). 



26 



[10] R Development Core Team: R: A language and environment for statistical computing. 
R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0 (2007), 
http : //www . R-pro j ect . org. 



[11] See, e.g., sec. 8.3 in R. Crandall and C. Pomerance: Prime Numbers: A Computational 
Perspective. Springer- Verlag New York (2001). 

[12] C. Asensio and A. Segui, in preparation. 

[13] D. Schwartz-Perlov and A. Vilenkin: Probabilities in the Bousso-Polchinski multiverse. 
JCAP 0606, 010 (2006), hep-th/0601162. 



27 



